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A THEOREM  ON  OPTIMAL  LOCATION 


A.J.  Goldiiian-  and  P . R . Meyers 
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ABSTRACT 


The  problem  studied  concerns  the  location  of  a facility  (e.g.  a 
mail-sorting  installation  or  a transport  terminal)  so  as  to  minimize 
the  sum  of  transport  costs  for  the  "source" -to-facility  and  facility - 
to-"sink"  flows.  Movements  are  assumed  to  occur  along  least  costly 
paths  in  a network  of  roads  or  other  conduits.  The  only  assumption  on 
cost  structure  is  that  marginal  cost  is  a non-increasing  function  of  trip 
length,  this  function  possibly  changing  from  trip  segment  to  trip  segment. 
The  distributions  of  sources  and  sinks  along  the  roads  can  be  arbitrary, 
as  can  the  network  topology.  Under  further  mild  hypotheses,  it  is  shown 
that  the  search  for  an  optimal  facility  location  can  be  confined  to 
endpoints  of  roads  and  interchanges  between  roads.  A detailed  study  is 
made  of  a radially  symmetric  circular-disk  model  in  which  only  radial  and 
circumferential  movement  is  permitted  (each  with  a constant  unit  cost), 
and  the  "density"  of  sources  and  sinks  drops  exponentially  as  a function 


of  distance  from  the  center  of  the  disk. 
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A THEOREM  ON  OPTIMAL  LOCATION 


(1) 


by 


A.J.  Goldman  and  P.R.  Meyers 
National  Bureau  of  Standards 

1.  INTRODUCTION 


Suppose  given  a collection  of  "customers"  for  the  services  of  a central 
facility(mail  sorting  installation,  transport  terminal,  firehouse,  library, 
school,  food  distribution  point,  etc.).  The  common  feature  of  the  class  of 
applications  we  have  in  mind  is  that  performance  of  the  service  involves 
the  movement  of  something  between  the  customer  location  and  the  central 

facility  e.g.  of  mail  from  patrons  to  sorting  installation  and  from 

sorting  installation  to  recipients,  or  of  travellers  from  initial  origins 
to  transport  terminal,  or  from  terminal  to  final  destinations. 

Thus  each  customer  location  can  act  either  as  a source  of  flow  to  the 
facility,  or  a sink  of  flow  from  the  facility,  or  both.  We  can  weight  each 
source  location  according  to  the  amount  of  flow  emanating  from  it,  and 
each  sink  location  according  to  how  much  flow  it  receives,  thus  giving  rise 
to  two  mass  distributions , denoted  jj,  (for  the  sources)  and  l#  (for  the  sinks) . 

The  flow  is  required  to  move  along  a least  costly  among  the  "allowable 
paths";  the  latter  will  be  called  roads  for  concreteness,  though  interpretations 
as  pipelines,  pneumatic  tubes  or  other  types  of  conduit  are  equally  in  order. 

We  use  the  notation  equally  in  order.  h 

( 2 ) 

c(P,L)  = cost  of  transport ingv  1 one  unit  from  source 
location  P to  facility  location  L , 

c*(L,P*)  = cost  of  transporting  one  unit  from  facility 
location  L to  sink  location  P*  . 

Then  the  total  transport  cost,  including  both  incoming  (to  the  facility)  and 
outgoing  flow,  is 

f(L)  = c(P^L)  jLt  (P)  + I c*(L,P*)  I#  (P*)  • (1) 

(1)  Supp or ted  by  the  Post  Office  Department  and  by  the  Department  of  Commerce's 
High  Speed  Ground  Transportation  Project.  No  official  endorsement  is 
implied. 

(2)  Transport  costs  are  considered  exclusive  of  the  costs  at  the  central 
facility  of  unloading,  processing  (servicing)  and  loading.  These  costs 
are  assumed  independent  of  the  location  of  the  facility. 


The  "optimal  location  problem"  for  a single  facility  can  be  defined,  for 
present  purposes,  as  that  of  choosing  the  facility's  location  L so  as  to 
minimize  f(L).  It  may  be  useful  to  "smear  out"  the  representation  of 
small  customers  to  a continuous  mass  distribution,  retaining  some  discrete 
mass-points  or  mass-curves  as  required  to  represent  extra- large  customers 
or  concentrations  of  them.  Thus  (l)  will  be  replaced  by 


In  recent  (unpublished)  work^  , C.  Witzgall  discussed  a number  of 
topics  relating  to  this  problem.  As  one  of  several  concrete  illustrative 
models,  he  considered  the  case  of  identical  uniform  distributions  of  sources 
and  sinks  over  a circular  disk.  Only  radial  and  circumferential  movements 
were  permitted,  a situation  suggestive  of  a city  which  has  a radial  web  of 
major  streets  and  is  surrounded  by  a beltway.  The  two  parameters  of  the 
model  were  both  costs  of  transporting  a unit  over  a unit  distance,,  namely 
the  costs  for  radial  and  for  peripheral  movements,  which  will  be  denoted 
c_  and  c^  respectively.  Witzgall  showed  that  the  condition  for  a peripheral 

location  of  the  facility  to  be  better  than  a central  one  is 


according  as  short-cuts  through  the  center  of  town  to  a peripheral  facility 
are  or  are  not  permitted. 

No  proof  was  given,  however,  that  the  optimum  had  to  occur  for  one  of 
these  extreme  cases  (location  at  center  or  on  circumference)  rather  than 
at  an  intermediate  position  along  a radius.  This  certainly  seemed  plausible 
(and  turned  out  to  be  true)  for  the  simple  disk  model,  but  its  main  interest 
lay  in  the  likelihood  that  it  represented  a special  case  of  some  general 
principle . 

The  present  paper  records  our  current  under standing  of  this  principle. 

As  will  be  seen,  it  admits  arbitrary  mass  distributions  a and  u*,  and  road 
networks  much  more  general  than  that  of  the  disk  model.  The  unit  (per- mile) 
transport  cost  along  a road  need  not  be  constant;  indeed,  the  critical  assum- 
ption describing  the  principle ' s range  of  validity  appears  to  be  that 
transport  cost  along  a road  (per  unit  carried)  is  a concave  function 
(perhaps  differing  from  road  to  road)  of  the  distance  Lvelled.  That  is, 


long-haul  costs  no  more  per  mile  than  u.oes  short-haul 

(3)  C.  Witzgall,  "Optimal  Location  of  a Central  Facility;  Mathematical  Models 
and  Concepts,"  1964. 

(4)  This  statement  about  average  costs  should  really  be  replaced  by  the 
corresponding  one  about  marginal  costs. 


f(L)  = r c(P,L)  dM  (p)  + J*  c*(L,  P*)  du*(p*)  . 


(2) 


c /c  < 0.216  , or  c /c  <2/3 it  = 0.212  , 


(3) 


2 - 


In  the  following  Section  2,  we  formulate  the  present  version  of  the 
principle  as  a theorem,  which  is  then  proved  in  Section  3.  Some  extensions 
are  described  in  Section  4.  Section  5 contains  an  analysis  of  a more 
general  "circular  disk"  model  with  radial  symmetry;  advantage  is  taken  of 
the  theorem,  in  that  the  search  for  an  optimal  location  for  the  central 
facility  is  confined  to  the  center  and  periphery  of  the  disk. 
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2.  FORMULATION  OF  THEOREM 


We  consider  two  systems  F and  p*  of  roads;  the  first  carries  the- 
source-to-facility  flow,  and  the  second  the  flow  from  facility  to  sinks. 

Some  roads,  or  segments  of  roads,  may  be  common  to  p and  P*  . The 
set  of  points  covered  by  the  roads  of  P will  for  simplicity  be  identified 
with  the  set  of  possible  locations  of  sources'''5'*  (the  distinction  between 
U > 0 and  t2=0  implicitly  sorting  out  the  "actual"  locations  among  the 
possible  ones),  and  this  point- set  is  assumed  compact (®);  similarly  for 
F*  and  the  possible  locations  of  sinks.  It  is  assumed  that  each  possible 
source  location  is  connected  within  P to  each  possible  location  of  the 
facility,  and  similarly  for  the  facility-to-sink  flow. 

Each  road  is  taken  to  be  a continuous  rectifiable  curve  in  the  usual 
mathematical  sense;  this  implies  an  orientation''''*,  and  a running  "path 
length  parameter"  which  increases  from  zero  at  the  beginning  of  the  road 
to  some  finite  maximum  value  (the  length  of  the  road)  at  its  end.  Some 
roads  may  be  closed  curves,  i.e.  their  beginning  point  and  ending  point 
coincide  (as  for  a beltway) . 

As  noted  earlier,  transport  cost  along  any  one  road  (per  unit  carried) 
is  assumed  a concave  function  of  the  road-length  traversed,  this  function 
possibly  differing  from  road  to  road.  While  this  seems  plausible  when  expressed 
as  a comparison  of  long-haul  per  mile  transport  cost  with  that  for  short- 
haul  (the  former  not  exceeding  the  latter),  it  may  introduce  a bias  into 
the  choice  between  central  or  peripheral  location  for  the  facility.  Flow 
from  an  outlying  source  to  a central  location  may  in  its  progress  encounter 
successively  more  congested  zones,  and  therefore  succesively  higher  per 
mile  costs;  our  concavity  hypothesis  in  effect  assumes  that  this  effect  is 
at  least  counter-balanced  by  the  economies  of  long-haul  operations.  On  the 
other  hand,  the  reverse  holds  true  for  flow  from  central  facility  to  outlying 
sink.  Roughly  speaking,  the  bias  is  likely  to  be  for  or  against  a central 
location  according  as  j u.  or  u*  is  more  concentrated  about  the  "center" 
of  the  region  under  study. 


(5)  This  evades  any  need  to  provide  a more  or  less  arbitrary  description  of 
how  flow  from  off-F  sources  moves  toward  p.  Such  a description  can  be 
guided  by  reality  for  any  specific  application. 

(6)  The  desire  to  encompass  idealized  models  with  infinitely  many  customers 
and  roads  and  special  points  requires  more  technical  mathematics  than 
would  be  needed  for  discrete  models.  Readers  unfamiliar  with  these 
concepts  (e.g.  compactness)  can  ignore  them  without  serious  loss. 

A bidirectional  road  can  be  represented  as  two  oppositely  oriented  roads. 
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(7) 


Certain  special  points  on  the  roads  of  the  two  systems  are  distinguished. 
These  are  required  to  include  the  beginning  and  ending  points  of  all  roads 
which  are  not  closed  curves,  and  of  all  segments  shared  by  any  specific  set 
of  roads.  (Thus  the  "focus"  of  a radial  web  of  roads  would  be  a special  point.) 
They  are  also  required  to  include  all  interchanges,  i.e.  points  where  two 
roads  intersect  and  transition  from  one  to  the  other  is  permitted.  (For  the 
the  simple  disk  model  described  earlier,  this  means  that  all  points  on  the 
beltway  are  special.)  The  set  of  special  points  is  assumed  compact. 

From  this  last  assumption  we  see  that  an  allowed  location  for  the 
central  facility  is  either  a special  point,  or  else'-0-'  lies  on  a road 
segment  free  of  special  points.  In  the  first  case  the  special  point  must 
be  an  interchange  from  a road  of  J*  to  one  of  P*,  while  in  the  second  case 
the  road  segment  must  be  common  to  a road  of  P and  one  of  P*  . The  set  of 
allowed  locations  L is  supposed  compact;  thus  f(L)  , if  continuous, 
does  actually  attain  a minimum  value.  The  continuity  of  f(L)  is  ensured 
by  requiring  c(P, L)  and  c*(L,  P*)  to  be  continuous. 

As  noted  earlier,  flow  is  assumed  to  move  along  a least  costly  of  the 
allowed  paths.  More  precisely,  if  c(P,rr, L)  is  the  transport  cost  for 
unit  flow  along  the  allowed  route  ^ in  P from  P to  L , then 

c(P,L)  = inf  c(P,rr,L)  ; (4) 

rr 

similarly  for  c*(L,  P*)  . 

The  remaining  assumption  is  that  all  flow,  whether  to  or  from  the 
facility,  must  be  routed  so  as  to  pass  through  at  least  one  special  point. 

This  suggests  viewing  the  special  points  as  seats  of  some  recording  or 
check-in  operation,  or  perhaps  of  consolidation  and  break-bulk  operations. 

The  hypothesis  is  automatically  satisfied  when  the  facility  is  itself  at  a 
special  point.  It  is  automatically  satisfied  for  flow  from  or  to  a 
customer  location  on  a road  not  containing  the  central  facility,  since  such 
flow  must  pass  through  at  least  one  interchange.  Thus  the  assumption  is 
restrictive  only  if  the  facility's  location  L is  on  a road  segment  free 
of  special  points  except  for  its  endpoints,  and  even  then  it  requires 
"unnatural"  routing  only  for  the  (presumably  small)  fraction  of  flow 
originating  or  terminating  on  that  particular  segment. 


(8)  The  point  is  that  an  allowed  location,  if  not  itself  a special  point, 
cannot  be  a limit  of  special  points. 


THEOREM.  Under  the  above  hypotheses,  at  least  one 
special  point  is  an  optimal  location. 

The  important  implication  is  that  when  seeking  to  determine  an  optimal 
location,  one  can  simplify  the  task  by  restricting  attention  to  special 
points.  This  may  be  useful  not  only  in  placing  new  facilities,  but  also  in 
explaining  historic  patterns  of  location  of  "central  places".  In  connection 
with  this  last  remark,  it  should  be  noted  that  cities  can  be  regarded  (in 
part)  as  "service  facilities"  for  the  "customers"  of  their  surrounding  areas. 


3-  PROOF  OF  THEOREM 


As  observed  above,  the  hypotheses  of  the  theorem  ensure  the  existence 
of  at  least  one  optimal  location  Lq  . We  assume  that  Lq  is  an  ordinary 

(i.e.,  non- special)  point,  and  prove  the  theorem  by  showing  that  this 
implies  the  existence  of  an  optimal  special  point. 

Because  Lq  is  ordinary  and  the  set  of  special  points  is  compact,  L_ 

must  lie  on  some  road  segment  rr  whose  interior  contains  only  ordinary 

points,  and  which  lies  both  on  at  least  one  road  of  T and  on  at  least  one 
road  of  r*  • Consider  the  continuation  of  through  Lq  in  one  of 

the  two  possible  directions.  A "branch- point"  may  be  reached  at  which  rr 
loses  its  identity  because  some  road  containing  it  diverges,  or  some 
new  road  portion  merges  with  it;  such  a point  is  a special  point.  The 
beginning  or  end  of  a road  may  be  encountered,  or  an  interchange  may  be 
reached;  again  a special  point  must  be  present.  If  none  of  these  occurs, 
then  we  must  eventually  return  to  ; in  this  case  rr0  is  a closed 

curve  consisting  of  ordinary  points  only.  Since  this  curve  has  no  inter- 
changes, and  all  customer  locations  must  be  connectable  with  Lq  , it 

follows  that  all  possible  customer  locations  lie  on  rr  , i.e.  the  points 

of  all  roads  lie  on  nrQ  • This  is  impossible  because  all  flow  must  pass 

through  a special  point,  whereas  rr ^ lacks  special  points.  The  conclusion 

is  that  the  maximal  continuation  of  rrQ  through  in  either  direction 

ends  in  a special  point. 


From  now  on  we  reserve  the  symbol 
ordinary  points  containing  L , including  its  (special)  end-points  E 


crQ  for  the  maximal  segment  of 


and  E-^  . (The  possibility  of  a closed  curve  with  Eq 

special  point  offers  no  difficulty.)  The  points  of  rrQ 

L(l),  where  \ is  the  arc-length  parameter  along  n 

i.e.  L(0)  = E and  L(\  ) = E,  , where 

o v max  1 


E, 


as  its  sole 


will  be  denoted 
measured  from  E 

o 


X = length  of 
max 


7 - 


o 


It  will  be  shown  that  either  EQ  or  E±  is  optimal.  This  will  be  done 
by  proving  that  the  function 


g(\)  = f(L(\))  , (5) 


defined  for  0 < A,  < X , attains  its  minimum  at  an  endpoint  of  this 
— — max 

interval.  Since  L lies  on  rr  , we  have 

o o’ 


min  f(L)  < min  g(l)  = min  f(L)  < f(L  ) = min  f(L) 
L Lfa  _ ° L 


from  which  the  optimality  of  Eq  or  E^  will  follow. 


Let  Q be  the  set  of  all  special  points  on  all  roads  of  T which 

contain  rr  , and  assume  o < 1 1 . Consider  any  possible  location  P 

o max 

If  the  flow  from  P to  L(A,)  would  move  wholly  along  a single  road  of  £ , 
then  the  last  special  point  it  would  encounter  is  either  Eq  or  E^  , 
which  both  lie  in  Q . If  not,  the  flow  enters  the  last  road  it  traverses 
at  some  interchange  which  is  a point  of  Q . So  eq(4)  becomes 


c(P,L(A))  = inf  (c(P,  q)  + c(q,L(l))}  . 

qyQ 


(6) 


We  will  use  the  abbreviations 
c (a.)  = c( q, l( i) ) , 

Si 

cU,P)  = c(P,q)  + c (A.)  . 

Si  Si. 

The  appearance  of  c (A.)  in  (6)  refers  to  movement  along  a single  road, 
and  so  by  hypothesis  this  function  is  concave  in  \ . That  is,  for 

each  triple  (A.  , l+,t)  with 


0 < A. 


< 1 


< 1 
max 


9 
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o < t < 1 


(7) 


the  equality 


c^(tA.  + (l-t)A,+  ) > tc^(A,  ) + (l-t)c^(A,+)  (8) 

is  valid.  Since  c(P,  q)  is  independent  of  A,  , it  follows  that  for 
each  q^Q  and  each  P , c^(A.,  P)  is  concave  in  A. 

(9) 

This  will  be  used  to  prove  that  for  each  P , 
c(P,L(\))  = inf  {c  (A.,P):  qcQj  (9) 

Si 

is  concave  in  X . Let  (A,  , A,+,t)  be  as  in  (7).  For  any  > 0 , 
there  is  a q*-Q  such  that 

c(P,L(tA.”  +(l-t)l+))  + > c (tX~  + (l-t)\+,P)  . (10) 

Q. 

From  (8)  and  (9)>  we  have 

c (tA.“  + (l-t)l+,P)  > tc(  P,  L(  A,” ) ) + (1-t)  c(P,  L(A.+) ) . (11) 

Q. 

Combining  (10)  and  (ll),  and  recalling  that  6 was  arbitrary,  we  have 

c(P,L(t\“  +(l-t)l+))  > tc(P,  L(  A,” ) ) + (l-t)  c(  P,  L(  A,+  ) ) , (12) 

and  so  (by  definition)  c(P,  L(A,))  is  concave  in  A. 


(9)  For  completeness,  the  proofs  of  known  properties  of  concave  functions 
are  incorporated. 
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Inequality  (12)  can  be  rewritten  as 

c(P,L(t\"  + (l-t)k+))  - tc(  P,  L(  X,”  ) ) - (1-t)  c(P,L(\+))  > 0 . (13) 

A non-negative  integrand  yields  a non-negative  integral  with  respect  to 
any  mass  distribution  ju  ; applying  this  to  (13),  we  obtain 

J'c(P,L(tk"  + (l-t)  A.+)d/Li(  P)  > tJ'c(P,L(A,"))dy(P) 

+ (l-t)fc(P,L(\+))du(P)  . 

This  is  simply  the  assertion  that  the  function 

gl(\)  = Jc(P,L(X.))da(P) 

is  concave  in  A,  . By  eqs(2)  and  (5)>  g(A.) 
similar  term  for  the  facility- to- sinks  flow, 
concave  by  similar  arguments,  and  so  g(X.)  is 
(X  , A,+,t)  as  in  (7), 

g(t  A.  +(l-t ) A,+  ) > tg ( A,"  ) + (l-t  )g(  A,+  ) 

> t min(g(  A.  ),  g(A,+  ))  + (l-t)  min(g(  A.”  >,g(b+) ) 

- min  (g(xT),  g(k+))  . 

That  i s 

g(x)  > min  (g(A.~),  g ( A,+ ) ) 

By  continuity  this  remains  true 
as  desired, 

min  (g(X)  : 0 < A,  < 
from  which  it  follows  as  shown  above  that  Eq  or  E^  is  an  optimal  location. 


for  A,  < X < X 


for  A,  =0  and  X '=  X so  that, 

max 


= min  6Uk)I  , 


is  the  sum  of  g^(A,)  and  a 
The  latter  can  be  proved 
itself  concave.  Thus,  for 
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b . SOME  EXTENSIONS 


One  extension  of  the  theorem  arises  if  its  hypotheses  are  weakened, 
to  require  concavity  of  transport  costs  along  only  some  roads  and  road 
segments.  The  preceding  proof  can  easily  be  adapted  to  establish  the 
existence  of  at  least  one  optimal  location  which  is  not  an  ordinary 
point  on  any  of  those  roads  and  segments  for  which  concavity  is  postulated. 

The  description  of  a second  extension  requires  another  definition. 

A function  h(l),  defined  (say)  on  an  interval  o < X < X , is 

^ — ■ max 

called  strictly  concave  in  X if  for  all  (x~ ,X  ,t)  obeying  (7),  the 
strict  inequality 

h(t\“  + (l-t)x+)  > th(xT)  + ( 1- 1 ) h ( X+ ) 


holds . 

If  g(l)  is  strictly  concave  in  X , the  last  part  of  the  preceding 

proof  can  be  modified  to  show  that  the  minimum  of  g(\)  occurs  only  for 

X = 0 or  X = X . It  would  follow  that  only  special  points  can  be 

max  ■ 

optimal . 

One  might  hope  that  the  strict  concavity  of  g(l)  could  be  inferred 
if  the  hypothesis  of  concavity  of  transport  costs  along  any  road  were 
strengthened  to  strict  concavity.  This  will  not  always  work.  The  functions 
c (X, P)  will  each  be  strictly  concave,  but  the  infimum  of  an  infinite 
collection  of  strictly  concave  functions  need  not  be  strictly  concave 
(it  can  for  example  be  linear),  so  that  "strictness"  can  be  lost  in  the 
transition  to  the  function  c(P,L(A,))  defined  in  (9)  • 

It  is  however  true  that  the  minimum  of  a finite  collection  of  strictly 
concave  functions  is  strictly  concave.  If  a is  discrete,  then  the  passage 
from  c(P,L(x))  to  g^(l)  also  preserves  strict  concavity.  So  for  discrete 
models,  only  special  points  can  be  optimal  if  transport  costs  are  strictly 


concave . 


11  - 


5 . APPLICATION  TO  A CIRCULAR  DISC  MODEL 


In  this  section  we  illustrate  the  previous  material  by  examining 
in  some  detail  a generalization  of  the  circular  disc  model  mentioned 
in  Section  1.  The  mass  distributions  jLt  (for  sources)  and  M*  ( for 
sinks)  are  assumed  to  coincide,  as  are  the  road  systems  L and  D* 
and  the  associated  transport  costs.  Thus,  for  any  location  of  the 
facility,  the  cost  of  the  sources-to-facility  flow  will  be  exactly 
half  that  of  the  total  (two-way)  flow. 

Polar  coordinates  (r,0),  with  origin  at  the  center  of  the  disc, 
will  be  employed.  The  disc  itself  is  characterized  by  0 < r < 1 , 
i.e.  its  radius  is  (for  simplicity)  chosen  as  the  unit  of  length.  The 
mass  density  is  assumed  to  fall  off  e xp one nt i ally  from  the  disc-center 
in  a radially  symmetric  way  , i.e. 

d/i(r,d)  = K exp  (-ar)  r dr  dd 

where  K is  a normalization  constant,  r dr  d0  is  the  element,  of  area  in 
polar  coordinates,  and  a > 0 is  a parameter  representing  the  degree  of 
central  concentration  of  the  distribution.  The  limiting  case  a=0  corres- 
ponds to  the  uniform  distribution  studied  by  Witzgall  (op  cit ) , and  our 
study  generalizes  his  precisely  in  permitting  a > 0 . 

The  system  of  available  roads  is  taken  to  consist  of  the  periphery 
of  the  disk  and  (initially)  all  of  its  radii,  regarded  as  bidirectional. 

(lO)  This  assumption  is  a fairly  common  one.  See  for  example  H.K.  Weiss, 
The  Distribution  of  Urban  Population  and  an  Application  to  a 
Servicing  Problem,  Operations  Research  9 (l96l) , PP-  860-87^- 
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The  intersection  of  each  radial  road  with  the  periphery  is  taken  as  a 
transfer  point  and  also  as  an  end-point  of  that  road,  the  other  end-point 
being  the  center  of  the  disk.  Unit  costs  of  peripheral  and  radial  travel 
are  assumed  to  be  constants  (c  and  c^  , respectively);  thus  the  transport 
cost  along  any  road  is  a linear  (hence,  concave)  function  of  distance 
travelled.  It  follows  that  our  general  principle  indeed  applies,  to 
show  that  the  optimal  location  for  the  facility  is  either  at  the  center  of 
the  disc  or  on  the  periphery.  By  symmetry,  we  need  consider  only  one 
among  the  peripheral  locations,  and  we  choose  the  one  given  by  0=0  . 

First  suppose  the  facility  is  located  at  the  center  of  the  disc. 

Then  the  least  costly  path  from  any  point  to  the  facility  consists  of 
a single  radial  segment,  and  so  the  total  cost  (for  the  two-way  flow)  is 

fc  = 2 J1  r crrdM(r,0) 

O -77 

= b 77  Kc  r^  exp  (-ar)dr 

o 

__  O p 

= U 77  Kc^a  “>[2-(2+2a+a  )exp(-a)]  . 

Second,  suppose  the  facility  is  located  on  the  periphery,  at  the 
point  given  by  0=0  . There  are  really  two  cases  to  consider,  depending 
on  whether  or  not  short-cuts  through  the  central  district  are  permitted, 
i.e.  whether  or  not  the  center  of  the  disc  is  assumed  to  be  an  interchange 
point  between  distinct  radial  roads.  We  assume  for  the  moment  that  such 
short-cuts  are  forbidden,  so  that  the  least  costly  path  from  any  point 
consists  of  a radial  trip  to  the  periphery,  followed  by  a trip  to  the 
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facility  along  a minor  arc  of  the  periphery.  The  corresponding  total  cost  is 


f = 2 J1  J77  [c  (1-r)  + c 1 0 1 ] d/i(r,0) 

^ O -77  ^ 

= b Kc  (l-r)r  exp  (-ar)dr 

o 

+ ^-Kc  r exp  ( -ar)  d d0  dr 

P o o 

= U-TT  Kc  a ^[(a-2)  + (a+2)  exp  (-a)] 


O _p 

+ 2 IT  Kc^a  [l-(l+a)  exp  (-a)]  . 


Our  objective  is  to  determine  the  "critical  value"  T*  = T*(a)  , 
of  the  ratio 


T = C /c  , 
p r 


(lM 


which  divides  the  cases  for  which  the  peripheral  location  of  the  facility 
is  optimal  from  those  for  which  the  central  location  is  optimal,  i.e. 

f < f if  T < T*( a)  , 

pc  ’ 

f > f if  T > T*(a)  . 

pc 

Thus  T*  can  be  found  by  equating  f to  fc  , and  solving  for  T . 

The  result  (for  the  case  of  forbidden  short-cuts)  is 

T*(a)  = (2/ff) [(4-a)exp(a)  - (4+3a+a2) ]/a[exp(a)  - (l+a)  ] , (15) 
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which  for  small  a has  the  form 

r*(a)  = (2/77)  [l+o(a)  ]/[3+a+o(a)  ] . 

Thus  r*(o)  = 2/3'77  = 0.212+  , in  conformance  with  Witzgall's  result  for 
the  uniform  distribution. 

Increasing  values  of  a correspond  to  greater  central  concentration, 
and  therefore  to  a greater  tendency  for  the  central  location  to  be  optimal. 
We  should  therefore  expect  T*(a)  to  be  a decreasing  function.  This  can 
indeed  be  verified  analytically  (details  omitted),  and  is  also  indicated 
by  the  second  column  of  Table  1.  It  is  especially  interesting  to  note 
that  T*(a)  turns  negative  at  roughly  a=2.7,  representing  a degree  of 
concentration  at  (and  beyond)  which  subsidization  of  circumferential 
travel  would  be  required  to  overcome  the  cost  advantage  of  the  central 
location. 

Next  we  assume  that  short-cuts  through  the  center  are  permitted. 

It  is  easily  seen  that  the  least-cost  path  from  any  point  (r,d)  to  the 
facility  is  either  the  same  one  assumed  in  the  previous  discussion,  or 
the  one  involving  a radial  trip  to  the  center  followed  by  a radial  trip 
along  the  ray  6=0  to  the  periphery.  The  former  is  less  costly  than  the 
latter  if  and  only  if 

c (l-r)  + c |6|  < c (l+r)  , 
r P r 

i.e.  if  and  only  if 


- 15  - 


' - 


TABLE  1:  Values^ 

of  t*( a) 

T*(a) 

A.  Short  cuts  through 

B.  Short  cuts  through 

a 

center  forbidden 

center  permitted 

0.05 

.208 

1 

.214 

0.1 

.205 

.210 

0.2 

.198 

.202 

0.3 

.191 

.194 

0.4 

00 

00 

1 — 1 

.186 

0.5 

.176 

.179 

0.6 

.168 

.171 

0.7 

.161 

.163 

0.8 

•153 

•155 

0.9 

IPs 

1 — 1 

.147 

1.0 

.137 

I 

.139 

1.2 

.121 

.123 

l.h 

.105 

.106 

1.6 

o\ 

co 

0 

.090 

1.8 

.073 

.073 

2.0 

.056 

•057 

2.3 

.032 

.032 

2.6 

.007 

.007 

2.7 

-.001 

.000 

(*)  Recall  that  T*( a)  is  the  value  of  c / c for  vhich  the  transport  costs 
associated  with  peripheral  and  central  facility  locations  are  equal. 


- 15a  - 


where 


2 2 

F(x,a)  = ax  + x(4+3a+a  )exp(-a)  - x'~(l+a)exp(  -a) 
- [6-6exp(-x)  - 2xexp(-x)]  . 


(18) 


Before  considering  this  equation  in  detail,  we  will  dispose  of  the 
remaining  case  T >2 flT  , corresponding  to  x > a . In  this  case  the 
short-cut  fails  to  be  best  if  and  only  if 

T 1 0 1 /2  < r < 1 , 

which  can  only  occur  for  | 0 | < 2/t  . Thus  the  total  cost  is 
f = 2 j‘2//7"  JT9//2  c (l+r)dju(r,e) 


o 

O 

+2  f/T 

J1  [c  (l-r)+c  d]d/i(r 

o 

70/2  r P 

+2  J 

J1  c (l+r)d^(r,0) 

2/T 

° r 

o 

: 2 Kc 

T 

(if57  j ^ r exp(-ar)dr  d0 

o o 

* f/T . 

j T^/2  r2  exp(-ar)dr  d0 

o 

o 

- J2/T  . 

~|  ^ 

(r0 r-r  )exp(-ar)dr 

o 

70/ 2 

+ r , 

r2  exp(-ar)dr  dd]  . 

2/T 

o 
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We  can  also  write 


f = 2 Kc  {J2/T  JT0/2  r2  exp( -ar)dr 


d0 


o o 


+ J2/T  r2  exp(-ar)dr  d0 

o 70/ 2 


J77  J1  r2  exp(-ar)dr  d0}  , 
2/t  o 


and  so 


(f  -f  )/2Kc  = j1  r exp( -ar)dr  d0 

^ o o 


+ J2/T  j2  (r0r-2r2)exp( -ar)dr  d0 
o 70/2 


- J"'  J77  r exfj( -ar)<10  di 


+ (70r-2rc  )exp( -ar)d0  dr 

o o 

= J1  [tt  + fr/T  (70-2r)d0j  r exp(-ar)dr 

o o 

= J2  {tt  - 2r2/7]  r exp(-ar)dr  , 

o 


which  has  a positive  integrand.  It  follows  that  for  7 > 2 /iT  , a 
central  location  is  always  preferable. 
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Thus  it  suffices  to  consider  (17)  and  (l8),  for  x < a . Here  (17) 

has  an  extraneous  root  x=0  for  all  a , arising  from  a multiplication  by  T 
in  simplifying  to  reach  (l6) . The  remaining  root  x*(a)  was  determined 
numerically,  and  then 

T*(  a)  = 2x*(a)/77a 

was  calculated^"'-)  . It  appears  as  the  third  column  in  Table  1.  We 
see  that  the  availability  of  the  short-cut  has  no  great  effect  on  T*( a)  ; 
it  is  again  a decreasing  function  which  turns  negative  at  approximately 

a=2 .9 • 

So  far  our  illustrations  have  dealt  exclusively  with  a continuous 
mass  distribution.  To  help  assess  the  consequences  of  such  an  idealization, 
we  conclude  by  considering  a more  discrete  version  of  the  circular  model. 

Here  the  roads  consist  of  N radial  "spokes"  uniformly  spread  around  the 
disk,  lying  along  the  rays 

d = k(27r/N)  , k = 0,1,  . . . ,N-1  , 

and  also  the  circumference  itself.  Along  each  radius,  the  mass  distribution 
(as  before)  decays  exponentially  as  a function  of  r.  For  simplicity  we 
treat  only  the  case  in  which  short-cuts  through  the  center  are  forbidden. 

By  our  general  theorem,  the  optimal  location  for  the  facility  will 
lie  either  at  the  center  of  the  disk,  or  at  the  tip  of  one  of  the 

(11)  The  assistance  of  Miss  S.  Rose  (NBS  Technical  Analysis  Division) 

in  the  algebraic  and  numerical  calculations  is  gratefully  acknowledged. 
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N spokes  (which  we  take  to  he  the  one  at  0=0.)  For  a central 


facility,  the  total  cost  of  the  flow  to  and  from  it  is 


_-i  p 

f = 2NKc  Jr  exp(  -ar)dr  , 
o 

which  is  n/277  times  the  corresponding  expression  for  the  continuous  model. 


For  a peripheral  facility,  the  cost  if  N is  even  is 


f = 2NKc  (l-r)r  exp(  -ar)dr 

P r o 

(H-2)/2 

+ 2Kc  (27r/N)(N/2+2  L k)  J r exp(-ar)dr 
P k=0  o 


= 2NKc  J"1  (r-r^)  exp(-ar)dr 

o 


+ 77HKc  r exp(-ar)dr  , 

P o 


which  is  again  N/2ZT  times  the  corresponding  expression  for  the  continuous 
model.  This  difference  from  the  continuous  model  is  of  course  simply 
due  to  failure  to  "rescale"  by  changing  K ; the  conclusion  (for  even  N) 
is  that  T*( a)  is  exactly  the  same  as  for  the  continuous  model. 


When  N is  odd,  the  cost  for  a peripheral  facility  becomes 


f = 2NKc  j'*-  (l-r)  r exp(-ar)dr 
P r o 

(K-!)/2 

+ 2Kc  (277/n)2  L k J r exp(-ar)dr 
p k=0  o 
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+ 


2NKcr  f1  (r-r2)exp(  -ar)dr 

o 


7r(N-N  ^)Kc  r exp(-ar)dr 

P o 


Thus 


(f  -f  )/2MKc  = J1  (r-2r")exp(-ar)dr 
p o 


+ (7TT/2)(1-N~2)  J1  r exp(  -ar)dr  . 

o 

For  the  continuous  case,  on  the  other  hand,  we  find  from  previous  formulas 
that 

(f  -f  ) /UlTKc  = (r-2r2)exp( -ar)dr 

pc  r 0 

+ (>n-/2)  J r exp(-ar)dr  . 

o 

-2  -1 

It  follows  that  T*( a),  for  the  present  model,  is  (l-N  ) times  its 
value  for  the  continuous  model. 

These  results  give  at  least  some  reassurance,  that  using  a 
continuous  idealization  in  such  problems  leads  neither  to  qualitative  nor 
to  unacceptable  quantitative  errors. 
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